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CONVERSION FACTORS, U.S. CUSTOMARY TO METRIC (SI) UNITS OF MEASUREMENT 


U.S. customary units of measurement used in this report can be converted to 
metric (SI) units as follows: 


Multiply by To obtain 
inches 25.4 millimeters 
2.54 centimeters 
square inches 6.452 square centimeters 
cubic inches 16.39 cubic centimeters 
feet 30.48 centimeters 
0.3048 meters 
square feet 0.0929 square meters 
cubic feet 0.0283 cubic meters 
yards 0.9144 meters 
square yards 0.836 square meters 
cubic yards 0.7646 cubic meters 
miles 1.6093 kilometers 
square miles 259.0 hectares 
knots 1.852 kilometers per hour 
acres 0.4047 hectares 
foot-pounds 1.3558 newton meters 
millibars 1.0197 x 10°93 kilograms per square centimeter 
ounces 28.35 grams 
pounds 453.6 grams 
0.4536 kilograms 
ton, long 1.0160 metric tons 
ton, short 0.9072 metric tons 
degrees (angle) 0.01745 radians 
Fahrenheit degrees 5/9 Celsius degrees or Kelvins! 


ITo obtain Celsius (C) temperature readings from Fahrenheit (F) readings, use 
formula: C = (5/9) (F -32). 


To obtain Kelvin (K) readings, use formula: K = (5/9) (F -32) + 273.15. 


SYMBOLS AND DEFINITIONS 


A elliptic function 

B elliptic function 

C wave phase speed 

E complete elliptic integral of the second kind 
H wave height 

He unrefracted deepwater wave height 
K complete elliptic integral of the first kind 
Ky, refraction coefficient 

K, shoaling coefficient 

L wavelength 

Le cnoidal wavelength parameter 

Ly linear wavelength 

Q littoral transport rate 

S beach slope 

T wave period 

U Ursell parameter = HL?/d? 

V(u,v) particle velocity 

d Stillwater depth 

fy cnoidal shoaling function 

g gravity 

k wave number 

p pressure 

t time 

X-y cartesian ordinates 

a angle of wave incidence 

n free-surface elevation 

p density 


o(x,z,t) velocity potential 


b,o subscripts which refer to breaking or deepwater values, respectively 


TRANSFORMATION OF MONOCHROMATIC WAVES FROM DEEP 
TO SHALLOW WATER 


by 
Bernard Le Mehaute and John D. Wang 


I. INTRODUCTION 


An understanding of many nearshore phenomena relies on the ability to 
predict the local wave climatology, given a deepwater wave description. 
For example, a quantitative description of longshore sediment transport is 
based on a knowledge of the wave characteristics in the surf zone. This 
report presents methods for determining the changes in the characteristics 
of a wave traveling over a variable bottom from deep water to shallow 
water. 


The acute sensitivity of the rate of littoral transport to wave 
breaking characteristics implies an accurate determination of these 
characteristics. The problem has numerous facets: 


(a) Given a deepwater unidirectional monochromatic wave, what 
are the breaking wave angle, depth of breaking, breaking wave height, 
and related quantities? 


(b) Given a multidirectional deepwater incident wave spectrum, 
what is the distribution of breaking wave characteristics and the 
"equivalent"' monochromatic wave used to determine the littoral drift? 


(c) How should a synoptic wave climatology be treated in order 
to determine the rate of littoral drift and related quantities? 


Only the first problem is addressed in this report. The relevant 
literature is reviewed, and a new hybrid wave theory is proposed to 
determine wave breaking characteristics on a sloped plane beach. 


II. NONLINEAR WAVE TRANSFORMATION 
il Nonlinear Wave Shoaling.- 


It is generally assumed that the wave motion over a gentle slope 
is the same as that on a horizontal bottom, and that there is no re- 
flection nor wave profile deformation. The wave motion is then 
determined so that the rate of transmission of energy or energy flux 
is constant over varying depth. 


The average energy flux through a vertical plane of unit width 
perpendicular to the wave propagation is 


5 t+T n p LD 
Bee =e J j } (gz + ea =x V ) udzdt (i) 


where 


oe = density 
t = time 
T = wave period 


d = water depth 


n = free-surface elevation 
g = gravity acceleration 
p = pressure 


V(u,v) = particle velocity 
Z = vertical ordinate 


In the general case, linear or nonlinear, where the flow motion 
can be expressed by a potential function $(x,z,t), the Bernoulli 
equation yields 
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in which case ¢ can be expressed at any order of approximation, such as 

given by a Stokesian power series. Even though classical solutions for 

cnoidal waves are irrotational, the potential function is not expressed 

but rather the solution for (n,u,v) is given; therefore, the energy flux 
for cnoidal wave is determined from equation (1) where (V“ = u“~ + w*). 


The results of all the calculations pertinent to linear wave 
theory and linear wave shoaling are given in Le Mehaute (1976). 


Instead of expressing > at a first order of approximation as in 
the linear wave theory, $ is expressed at a higher order in equation 
(3), the shoaling coefficient Ko = H/Hg becomes not only a function 
of d/L or d/Ly but also a function of the deepwater wave steepness, 


Bley ry 


This calculation has been performed at a third order of approxi- 
mation (Le Mehaute and Webb, 1964), and the fifth order of approximation 
(Koh and Le Mehaute, 1966) based on the third-order solution and fifth- 
order solution for a Stokesian wave as developed by Skjelbreia and 
Hendrickson (1960). The first definition of Stokes for the phase velocity 


is used; i.e., the average horizontal water particle velocity over a 
wavelength is zero. The results of such investigation are presented 
in Figures 1 and 2. 


The correction AH due to nonlinear effects never exceeds 5 
percent and is more commonly of the order of 1 percent. These in- 
vestigations show that: 


(a) The nonlinear shoaling coefficient is initially less than 
the linear coefficient when d/Lo > 0.4, then becomes larger toward 
shallow water until the wave breaks. 


(b) The Stokesian power series is not uniformly convergent, i.e., 
the function of d/L of higher order tends toward infinity when d/L tends 
to small values. Therefore, the "best" order of approximation is not 
necessarily the highest order. For relatively deep water d/L > 0.25, the 
fifth order of approximation would be the best insofar as wave height 
transformation is concerned; for very shallow water d/L < 0.01, the 
linear theory would be best. In the intermediate range the third-order 
theory would be best, and therefore should be preferred overall because 
of its range of applicability. 


The second definition of Stokes for the phase velocity can also be 
used; the average momentum over a wavelength is zero by addition of a 
uniform motion. Yamaguchi and Tsuchiya (1976) indicate that the results 
yield slightly larger values, at most a 7-percent increase for the 
shoaling coefficient, than the results obtained by Le Mehaute and Webb 
(1964). 


The principle of conservation of energy flux has also been applied 
to a cnoidal wave, and like the Stokesian wave the results depend on the 
order of approximation and the definition of phase velocity. All these 
investigations on cnoidal waves are based on an energy flux such as ex- 
pressed by equation (1). Masch (1964) was the first to deal with this 
subject; however, his wave theory is not consistent, even erroneous, 

(in the table of functions used by Masch in the shoaling of cnoidal 
wave, the water depth below MWL should be substituted by ht, the water 
depth under trough), and the results are presented in a form which is 
difficult to use. The relation to deepwater wave and sinusoidal theory 
is not discussed and no attempt is made to follow the shoaling of a 
specific wave. 


A significant contribution to the shoaling of cnoidal waves is 
given by Iwagaki (1968). Iwagaki treats the case of an approximate 
solution of cnoidal wave in which he used the second definition of 
phase velocity as given by Laitone (1961). The approximation is on the 
value of the elliptic integral which is replaced by a simple function 
of empirical coefficients. Iwagaki shows that this simplification 
actually covers a wide range of cases and allows him to simply in- 
vestigate the shoaling of what he calls "hyperbolic waves."' When the 
energy flux in deep water (as computed using small-amplitude theory) 
is equated to the energy flux in shallow water, described by first- 
order hyperbolic waves, Iwagaki obtains 
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Figure 1. The shoaling coefficient at the third order 
of approximation (Le Mehaute and Webb, 1964). 
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Figure 2. The shoaling coefficient at the fifth order 


of approximation (Koh and Le Mehaute, 1966). 
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According to Iwagaki (1968), this theory yields sufficiently accurate 
results for Ursell parameter U > 47. However, as pointed out by 
Svendsen (1974), the theory of Iwagaki deserves to be regarded as a 
practical solution to second-order cnoidal waves when the deepwater 
wave steepness is smaller than 0.02 and the relative water depth is 
smaller than 0.05. The matching of the Iwagaki hyperbolic wave with 
the third-order Stokesian wave is shown in Figure 3. 


The shoaling of the true cnoidal wave has been investigated by 
Svendsen and Brink-Kjaer (1972), Svendsen (1974), and Svendsen and 
Hansen (1977). They also give H/Hj as function of d/Lo and Ho/Lo 
(Fig. 4) and a computer-printed table. It can then be shown that for 
large values of Ursell parameters the shoaling coefficient K, + d-1 
instead of d-1/4 as given by the Green law (long wave linear theory). 
Concurrently, Shuto (1974) arrives at very similar results. 


Yamaguchi and Tsuchiya (1976) also carry out the same calculation 
based on the two definitions of the Stokes wave velocity for the cnoidal 
theory of Laitone (1961) and that of Chappelear (1962). However, an 
arithmetic error has been found in the Laitone theory (Le Mehaute, 1968). 


Zo Comparison and Matching Between Various Theories. 


As a wave propagates from deep water to shallow water it is. 
theoretically possible to determine the variation of wave height, wave- 
lengths, etc. This could be done by applying the principle of con- 
servation of energy flux to either the linear wave or the nonlinear 
Stokesian wave, or the cnoidal and solitary wave. Since a Stokesian 
wave rather applies in deep water, the transformation of water wave 
should be followed with that theory for the largest value of relative 
depth d/Lo and then switched to the cnoidal theory when d/Lo becomes 
small. However, such a scheme implies that the theories can be matched 
continuously, but there is a priori no reason why the ratio H/Hj should 
be the same for the value d/L, which corresponds to the limit of validity 
of both theories. On the other hand, if the wave heights are matched, 
then the energy flux will present a discontinuity (Fig. 5). The signifi- 
cant feature is that the cnoidal wave height grows faster with decreasing 
depth, although at intermediate depth its value is up to 10 percent less 
than predicted by a Stokesian theory. Waves with wave steepness larger 
than 2 to 3 percent will break at a depth where the cnoidal wave height 
is only slightly larger than that of a Stokesian wave. Waves with small 
wave steepness, however, such as swells, reach much smaller depth before 
they break and consequently a major part of their shoaling process is 
governed by the cnoidal wave theory. For these waves, the two theories 
such as the Stokesian (first order or linear theory) and cnoidal wave 
at a second order will yield significantly different results. 
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Figure 3. Shoaling coefficient obtained in equating the 
energy flux between a third-order Stokian wave 
in deep water and cnoidal (hyperbolic) wave in 
shallow water (Iwagaki, 1968). 
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Figure 4. Comparison between shoaling coefficients 
according to linear and cnoidal wave theories 


(Svendsen and Brink-Kjaer, 1972). 
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Figure 5. Matching Stokesian (first order) and cnoidal 
wave theory (Svendsen and Hansen, 1977). 


These results (Fig- 5) show that no continuous transition is pos- 
sible between the two theories. This means that it is not possible to find 
a value of the water depth, d, where the curves for the two theories fit 
smoothly together. If the Stokesian theory is used in deeper water and 
changed to a cnoidal theory when the wave enters shallow water, there 
will be a discontinuity in the variation of either wave height or wave- 
length, or both, depending on which water depth is chosen for the switch. 
Of course, the same will appear for all other quantities such as particle 
velocities, pressure, etc., and the rate of change of these. Svendsen 
(1974) shows that the limit of applicability of the cnoidal theory is 
d/Lo < 0.1193 when H is small. Koh and Le Mehaute (1966) also showed 
that the limit of applicability of the fifth-order Stokesian wave theory 
1S} GV/Jh 2a Weil) when H/Lo = 0.05 and d/Lo > 0.13 when H/Lo = 0.10 (see 
Rie, 2). 


There is a large difference between Stokesian and cnoidal wave 
between d/Lg equal 0.1 and 0.3. In this region no known wave theory 
fits very well. It could have been expected that a higher order 
Stokesian theory would be the answer, but the investigation by Koh 
and Le Mehaute (1966) shows that when d/L, decreases the fifth-order 


approximation represents an even worse approximation than the third 
order. Similarly, it is found that second-order cnoidal theory is 
worse than first-order cnoidal theory for large wave steepness. This 
is inherent to the point that both cnoidal and Stokesian power series 
expansion in terms of the small parameters h/d and H/L respectively 
are nonuniformly converging series since the functions of d/L attached 
to each power term blow up when d/L tends toward small values. 


It is interesting that Yamaguchi and Tsuchiya (1976) found that 
the shoaling coefficient given by Le Méhaute and Webb (1964) (first 
definition of Stoke's phase velocity) almost coincides with the 
shoaling coefficient obtained from cnoidal theory developed by 
Chappelear (1962) (second definition). 


Shuto (1974) attempted to make a synthesis of all these theories 
in a simple and practical form by empirically matching these solutions. 
Subsequently, he proposes the following law for practical purposes: 
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These equations seem to be the most realistic to remember from aulpler 
the theoretical approaches. In the range where both cnoidal and third- 
order Stokesian theory apply, the values of the shoaling coefficient 
are very close to each other as shown in Figure 6 (Flick, 1978). 
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Figure 6. Comparison between Stokesian third order and 
cnoidal shoaling coefficient with experiments 
(from Flick, 1978). 


Interestingly, the use of the linear wave theory to evaluate the 
value of the shoaling coefficient extends much beyond the formal 
validity of this infinitesimal wave theory. Similarly, the value of 
the shoaling coefficient given by the Stokesian wave theory extends 
into the area where the cnoidal theory fits best. This is due to the 
fact that, the shoaling coefficient being the ratio of wave height 
H/Hg only, the increase in free-surface elevation under the crest is 
partly balanced by the increase of free-surface elevation under the 
wave trough. However, that the linear wave theory applies for the 
shoaling coefficient does not mean that all wave characteristics 
(wavelength, velocity components, pressure, acceleration) follow the 
same principle; after the local wave height is obtained, all other 
wave characteristics are determined by the appropriate theory. 


Se Comparison Between Theory and Experiment - 


A relatively large number of experiments have attempted to 
verify shoaling laws; all have been conducted in laboratory wave flumes 
with waves generated by wave paddle. Most of these experiments suffer 
lack of accuracy because they were either done at too small a scale 
and were subsequently subjected to significant scale effects such as 
large viscous damping experiments (Iversen, 1951), or the wave paddle 
generated not only monochromatic waves but harmonic components 
(solitons) which introduced significant error and scattering 
(Eagleson, 1956; Iwagaki, 1968). 


There is actually considerable controversy whether waves of 
steady-state profile exist, as demonstrated by Dubreuil-Jacotin (1934). 
Theorists Benjamin and Feir (1967) and experimentalist Galvin (1970) 
postulate that the disintegration of finite amplitude monochromatic 
wave occurs in deep water even on horizontal bottom. There are as many 
theoreticians who assume that a steady-state profile does exist as there 
are experimentalists who do not notice the "creations" of solitons. 


Use of a formulation developed by Mei and Le Mehaute (1966), 
Peregrine (1967), and Madsen and Mei (1969) indicates that for a 
sufficiently abrupt change in water depth, both a solitary wave and a 
cnoidal wave disintegrate into multiple crests. These results have 
been obtained numerically and verified experimentally. However, over 
a relatively gentle beach, the wave period remains constant between 
deep water and shallow water and no disintegration takes place. Dis- 
integration takes place when the wave arrives on a reef. It seems 
natural to assume that the difference between these two observations 
is due to the difference in bottom slope. Benjamin and Feir (1967) 
show that waves are unstable if kd > 1.4; however, experiments by 
Flick (1978) indicate that kd can be much larger without evidence of 
wave disintegration or spectral smearing. 


It is commonly accepted that a monochromatic wave arriving on a 
rapid change of depth (in diffraction zone) gives rise to at least a 
doubling of crests. Such phenomenon is due to the nonlinear con- 
vective effects. Iwagaki and Sagai (1971) also investigated the 


deformation of long waves over a gentle slope using the nonlinear long 
wave theory and power series expansions. They found the shoaling 
coefficient to be a function of beach slope when S > .01. The steeper 
the slope the smaller the shoaling coefficient, a fact which can be 
attributed to partial reflection. In fact, due to friction effect, the 
ratio H/Hg for a given value of d/Lp decreases instead of increases 
(Sawaragi, Iwata, and Masayashi, 1976). 


The first reliable experiments were conducted by Brink-Kjaer and 
Jonsson (1973) and Flick (1978). Figures 7 and 8 show results for different 
values of H,/Lj. Flick separates the first, second, and third harmonics 
from his wave data and is subsequently able to give a reliable ex- 
perimental shoaling coefficient. Flick compares his results with 
Le Mehaute and Webb (1964) (third-order Stokesian) and also with the 
cnoidal solution of Svendsen and Brink-Kjaer (1972) in shallow water 
(see Fig. 6). 


The shoaling coefficient of a hyperbolic wave is also fairly well 
verified by Iwagaki (1968) who gives results very close to the two 
mentioned above. 


Svendsen and Hansen (1977) compared the shoaling of cnoidal wave 
with a set of careful experiments and claimed that other experimenters 
(Wiegel, 1950, Iversen, 1951; Eagleson, 1956) carried out their ex- 
periments on too steep a slope for the shoaling theory to be valid. 
Furthermore, they calculate the damping due to viscous friction, 
obviously important on a gentle slope. Svendsen and Hansen concluded 
that if the wave height at depth d/L, = 0.10 is matched between cnoidal 
and linear, rather than the energy flux, the cnoidal theory predicts 
the shoaling quite well, even close to breaking with small deepwater 
wave steepness Ho/Lg < 3 to 4 percent but not beyond. Consistently, 
with all theories, the wave just before breaking suddenly peaks up very 
rapidly (Le Mehaute, 1971). In this range of values, all shoaling 
theories (third Stokes, cnoidal and hyperbolic) tend to slightly under- 
estimate the value of the shoaling coefficient. Subsequently, the cal- 
culated breaking wave height tends to be underestimated. The linear 
wave theory underestimates the breaking wave height most significantly, 
sometimes by a factor of almost 2 (Fig. 9). 


It is pertinent to remember that (a) the shoaling coefficient 
given by the linear theory is valid beyond the limit generally con- 
sidered as valid for a linear theory, and (b) the shoaling coefficient 
given by third-order Stokesian wave is fairly well verified ex- 
perimentally and actually very close to the value given for the cnoidal 
wave, even though, as in the case of the linear wave, free-surface 
profile, pressure, velocity, and acceleration could be significantly 
different. 


In general, the linear theory can be applied throughout from deep 
water to shallow water and then the linear breaking wave height is multi- 
plied by a coefficient function of the beach slope (Koh and Le Mehaute, 
1966). After the wave height, H, is determined as a function of the 
deepwater wave height, H,, and wave period, T (or deepwater wavelength 
Lo), all other shallow-water characteristics (free-surface profile, 
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Comparison of experimental shoaling coefficients with Stokes 
third order (xp is the breaking location) (from Flick, 1978). 
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particle velocity, acceleration, and pressure) follow by application 
of one of the classical wave theories within the accuracy which is 
determined for the chosen theory. 


BREAKING INDEX. 


LINEAR 


EXPERIMENTAL 
SQ SHOALING 
La 


SS 


Figure 9. The linear theory underestimates the breaking wave height. 


A Nonlinear Wave Refraction, 


It has been shown previously how to determine the shoaling co- 
efficient, K, = H/Hj, when the wave arrives perpendicular to the 
bottom contour. This discussion deals with the refraction coefficient, 
K,. Refraction occurs when a wave arrives at an angle a with bottom 
contours; then H/H, = Kg Ky. For a straight parallel contour Snell's 
law becomes 


SGT constant (6) 
which applies whether the wavelength, L, is expressed by a linear theory 


or not. 
1/2 WZ 
se _ fcos 
0. (a) 


Then 
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which also applies for nonlinear as well as for linear theory. The 
subscript o refers to deepwater wave characteristics. 


In many cases, the refraction method provides a reasonably 
accurate measure of the changes waves undergo on approaching a coast. 
However, if the angle of a wave ray with the bottom contour is large 
(i.e., larger than 70°), minor error in the value of the incident 
angle leads to a large error in direction angle a in shallow water. 
Also, accuracy as far as height changes are concerned cannot be ex- 
pected where bottom slopes are steeper than 1/10. No strict limit 
has been set, but the accuracy of wave heights derived from orthogonals 
that bend sharply is questionable. In short, refraction coefficients 
which are quite different from unity, such as K, < 0.5 and K, > 1.5, 
must be doubted (Whalin, 1971). 


Nonlinear effects, having an effect on wavelength, phase and 
group velocity and energy flux, subsequently have an effect on wave 
refraction. This problem has been examined by Chu (1975) who used a 
mix Of three theories, i.e., the first-order cnoidal theory of Korteweg 
and DeVries (1895), the second-order hyperbolic wave of Iwagaki (1968), 
and the Stokes third-order wave as given by Le Mehaute and Webb (1964), 
which led to some inconsistencies in approximations. Skovgaard and 
Petersen (1977) used instead the first-order cnoidal theory of Svendsen 
(1974) and the stream function wave theories of Dean (1970). 


Theoretically, it is possible to express phase velocities as a 
function of the relative wave heights from nonlinear wave theories. 
For example, the deepwater wavelength at a third-order Stokesian 
approximation and the breaking wavelength by a cnoidal or hyperbolic 
wave theory can be conveniently expressed. However, it is interesting 
that due to deformation of wave profile on a sloped bottom, the simple 
linear theory has been verified (experimentally) quite well (Ippen, 
1966). Wavelengths given by linear and cnoidal theories are compared 
in Figure 10. Although the cnoidal theory predicts wave height well up 
to breaking, it overpredicts wavelengths significantly. Cnoidal theory, 
in fact, predicts an increase in wavelength for a decrease in depth 
when the relative height, H/d, is sufficiently large. This increase 
is not reflected by known data (Ippen, 1966) which are fitted quite 
well by linear theory (Fig. 11). 


III. BREAKING WAVE CHARACTERISTICS ON A SLOPED 
PLANE BEACH 


Ne Review of Previous Work. 


The determination of longshore currents and sediment transport 
depends crucially on the characteristics of the breaking wave field. 
The wave energy transport, or energy flux, is of particular importance 
such that accurate determination of wave height, wavelength, depth at 
breaking, and breaking wave angle becomes essential. 


This section deals with the practical aspects of determining 


the breaking wave characteristics when certain deepwater character- 
Usties are gaven. [he objective is) to derive and present cesullitis 
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— Linear Theory 


-—-- Cnoidal Theory 
Cnoidal Third-Order Stokes 


Figure 10. Wavelength as a function of depth pre- 
dicted by different wave theories. 
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Figure ll. Transformation of wavelengths (from Ippen, 
1966; used with the permission of McGraw- 
Hill Book Company). 
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consistent with present knowledge and in a readily usable form. 


The general problem would require the determination of the 
shoaling and refraction of a multidirectional wave spectrum from deep 
water over a randomly varying bottom topography until breaking occurs. 
Although such an analysis is possible, it is much too complicated and 
would have to be dealt with on a case by case basis using either 
manual or computer methods. 


A significant and useful simplification is achieved by assuming 
the bottom to be a uniformly sloped plane. This allows bottom vari- 
ations to be described by a single parameter, i.e., the bottom slope 
S. The refraction process is then described globally by Snell's law. 
This discussion deals only with monochromatic waves under the usual 
assumption that the wave period remains constant, and that friction 
and reflection are ignored. 


To obtain accurate predictions it is necessary to have a wave theory 
which is applicable up to the point of breaking. A lack of knowledge 
of the actual breaking process requires the use of an empirical 
breaking criterion to determine the point of breaking. A study by 
Le Mehaute and Koh (1967) evaluated the Stokes first-, third-, and 
fifth-order theories and compared the Miche (1944) breaking criterion 
with an empirically derived equation. One of these equations was 
derived by fitting a number of experimental data points covering the 
range 0.02 < S < 0.2 and 0.002 < Hj/Lo < 0.09. This equation ex- 
plicitly accounts for beach slope and is 


7 -1/4 
b Wa fo) 

mee 0.76 S a (8) 
fe) Ce) 


Since equation (8) is based on observed data it takes into account 
nonlinear effects such as wave height peak-up just before breaking. 
In applying this equation to waves arriving at an angle to the shore, 
Le Mehaute and Koh (1967) corrected the bottom slope for the angle of 
incidence; however, they neglected to replace the deepwater wave 
height with its unrefracted value. 


Subsequently, a new and easier approach to compute cnoidal 
waves was presented by Svendsen (1974). Brink-Kjaer and Jonsson 
(1973) showed that near breaking the water depth is usually so shallow 
that cnoidal theory applies. Indeed, it has been found that wave 
height is described well by cnoidal theory in the area close to and 
before breaking. 


In a recent report, Ostendorf and Madsen (1979) propose to use 
cnoidal and linear Stokes wave theories in their respective areas of 
applicability. A transition between the two theories which assumes 
continuous variation of energy flux and phase velocities is also 
presented. Ostendorf and Madsen further suggest the use of an em- 
pirical breaking criterion which is sensitive to bottom slope and 
depth-varying wave parameters, i.e., 
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To obtain the breaking wave characteristics, the two offshore 
parameters (sin a/C*, Cy) must be known where 


a = angle of incidence 


C = wave phase speed C* = C/(gT) 


g = gravity acceleration 
T = wave period 
Cr, at Vg ier n sin oon 
d = stillwater depth 
1 2kd 
Nay Cl eran 
ik Ss Qy/it, 


L = wavelength 


In the deepwater limit this implies that wave height, Hg, angle 
of incidence, a9, and wave period must be known independently. The 
solution requires an iteration process and nomographs are presented 
to facilitate the operations. 


The method for determining breaking wave characteristics sug- 
gested by Ostendorf and Madsen (1979) has been compared with experi- 
mental data (Kamphuis, 1963), and it was found that the predicted 
breaking wave angle is too large, especially for smaller wave steep- 
nesses (see Table 1). This is easily explained when considering the 
plot of wavelength transformation shown in Figure 10. Although 
cnoidal theory predicts wave height well up to breaking, it overpre- 
dicts wavelengths significantly. Cnoidal theory, in fact, predicts 
an increase in wavelength and therefore, also in wave angle when H/d 
is sufficiently large. This increase is, as previously mentioned, 
not reflected by known data (Ippen, 1966; Fig. 11), which are fitted 
quite well by linear theory. As another consequence, the wave break- 
ing criterion, again, a result which is difficult to defend. 


Dean (1974) determined wave breaking angles using his stream 
function theory, but with a slope-independent semiempirical breaking 
criterion. A comparison of his results with the experimental ob- 
servations of Kamphuis (1963) is also presented in Table 1. The 


BS 


predictions are consistently too high, especially for smaller wave 
steepness, where predicted and observed ap differ by a factor of 
approximately 2. 


Table 1. Comparison of measured and predicted 
breaking wave angles. ! 


Ostendorf 
and Madsen 
(1979) 


Dean 
(1974) 


1Beach silope, S = .1 


ie Solution Approach. 


To obtain reliable prediction of breaking wave characteristics, 
this study proposes to use cnoidal theory to describe the transformation 
of wave height while wavelength will be transformed using either linear 
wave theory or third-order Stokes theory. The cnoidal wavelength is 
then considered as an auxiliary parameter which cannot be identified as 
the physical wavelength. Linear wave theory is simpler to use; however, 
to retain some nonlinear effect in the transformation of wavelength the 
third-order Stokes theory is also included. The wavelength computed 
using the cnoidal third-order Stokes theory is shown in Figure 11. 


Due to the large wave heights near breaking, a higher order approxi- 
mation to cnoidal waves given by Iwagaki (1968) was considered. A more 
detailed discussion of this "hyperbolic" wave theory is given in the 
Appendix. Note that this higher order theory suffers from the same prob- 
lem of inhomogeneous convergence as, for example, plagues fifth-order 
Stokes waves (Le Mehaute and Koh, 1966), and therefore gives poorer 
results than the first-order cnoidal theory near breaking. 


The computation and shoaling of cnoidal waves have been given by 
Svendsen (1974) and Svendsen and Brink-Kjaer (1972). It is convenient 
to define the parameter 


HL2 
Up pag mcm (10) 
5 S 
d 
where 
H = wave height 
Le = cnoidal wavelength parameter 
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d stillwater depth 


K(m) complete elliptic integral of the first kind 


The dispersion relationship may be written as 


L 
eG fa H 


- 8 72 cs Ween eo 
The deepwater wavelength Lo or T and A = A(m) = = 1 aK 


The complete elliptic integral of the second kind is designated E. 


Invoking energy flux conservation between wave rays and using 
linear theory in deep water, the wave height transformation is given 


as 
MIS Kx Ls i 
Le ae € mo) S) f (12) 
H, \i6 le be H 


In equation (12) Bn = £ CU) = uo? : ae 


and 


B = B(m) ts E (6n" s Sn > 2 2 (Gn o 2) ‘) 


-(-9-f (13) 


Equations (12) and (13) define the shoaling of cnoidal waves and are 
used to determine the shoaling coefficient. 


To compute the wavelength and refraction the Stokes wave theory 
is used. Linear waves are described by the dispersion relationship 


L 
2 
pale tanh Zieh (14) 
IL, L 
(0) 1 


For the third-order Stokes approximation from Le Mehaute and Webb (1964) 
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= 22 
co : (22) 


and the refraction coefficient 


/cosa, 
KR : cosa (23) 


where a = angle of wave with shoreline. Also, when refraction is included 
the deepwater wave height in the above expressions should be replaced by 
its unrefracted value. 


Finally, the point of breaking is determined by a modified form of 
the empirical breaking criterion (eq. 8). To adapt the formula to waves 
approaching at an oblique angle, the bottom slope and deepwater wave 
height are replaced by S cosop and Ho (cosag/cosap)!/2, respectively. 
The breaking criterion is then obtained by 


me Ye 
we 0.76 S) woghl oh, (2) 7 (24) 
or 
ae 0075 ease!” pec ee K a (25) 
Hy r “bD use r 


Rewriting equation (25) using H,/H, = KsKy, where K, is the shoaling co- 
efficient yields 
H -1/4 


2 ys oO angle! Le O (26) 


Equation (26) gives the form of the breaking criterion used in this study. 
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An explicit solution for the breaking wave characteristics cannot 
be obtained from equations (11) to (23) and equation (26). A numerical 
solution is required and it becomes important to reduce the number of 
independent parameters as much as possible. By straightforward manipu- 
lation of the equations, only three independent parameters need to be 
known: deepwater wave steepness, H)/Lo; beach slope, S; and deepwater 
incident wave angle, ao. 


A computer algorithm is constructed to solve the problem. The 
basic process consists of the following steps: 


H 
(1) Give values for > S, % 
Lip 
Lio 
(2) Assume KK sA=1 and or L3pf = 0.4 
L3o0 
Hb 
3 1 fe 
(3) Guess a value for ae 
(4) Find 
EoD i \ oe Ho Hh, 
L. = ame K_K as 27 1 + ae 
O b Seen 5 
H H -l 
(5a) Determine = K. K x 2 (72) 
1b z o “\“lo 
H,\-1 H 
b b 
and k dad. = 2m ie ees 
ae 8 oS Lip 
ai 
: 6 3b 
(Sb) Determine i Kier Kennex =e 
3b Sunt Bion S36 
H\ -! H 
and kd) = 2n (2) ae 
3b 
Eis 
(6a) Find Tce tanh Kip qd 
lo 
1+C_A 
L 
(GD) Hint 22 6 sean ten 168 
Ib; 3b b 
30 1+ 03 


(7) Find o, from equation (22), using wavelength 
b q : 
ratio from step (6), and find K, from equation 
(16) 


a8) 


2, 
L_b Ib 
(8) Determine U = ao BS = 


lee b 


(9) Find m and K(m) satisfying equation (10) using 
polynomial approximation of K given in Abramovitz 
and Stegun (1964) 


(10) Determine E, A, B, f,, and then K_ from equation 
(12) H iS 


(11) Go back to step (4) until Ky remains constant 


(12) Compare the obtained K. value with value computed 
£rom equation, (26). it wdatferent.. somback tomstep 


(3) 
Se Results. 


The wave breaking angle, ap, as computed with linear, cnoidal- 
linear, and cnoidal third-order Stokes approaches, is compared with the 
experimental data of Kamphuis (1963) as shown in Le Mehaute and Koh (1967) 
and Figure 12. These are apparently the only data on ap and are ob- 
tained for a single bottom slope, S = 0.1. For Ho/Lo = 0.0175, linear 
theory gives the best fit to the data. On the other hand, cnoidal 
linear theory provides a better fit for the larger steepnesses, Ho/Lo = 
0.053 and 0.062. Unfortunately, it is not possible to draw any de- 
finite conclusions from the limited data. However, it appears that 
linear theory provides the best estimate of wavelength irrespective 
of wave steepness, as found by Eagleson and Dean (Ippen, 1966). Also, 
for large wave steepness, cnoidal theory predicts the wave height quite 
accurately up to the point of breaking. 


For a given beach slope the waves break at increasing relative 
depth ratios, d/L, as the deepwater steepness increases. For large 
enough Hj)/Ly, the cnoidal theory is no longer valid since it predicts 
a nonphysical complex wave height (Svendsen, 1974). The critical 
deepwater wave steepness for which the cnoidal theory ceases to be 
valid has been determined for five different bottom slopes (Table 2). 
This critical value is only weakly sensitive to the magnitude of ao. 
For Ho/Lg greater than the critical value, a different wave theory 
such as Dean's (1974) or third-order Stokes (Le Mehaute and Koh, 1967) 
must be used. 


In Figures 13 to 17, the variation of breaking wave angles 
versus deepwater wave angle with H)/Lo and S as parameters is de- 
picted as computed using the cnoidal-linear theory. Similar results 
for linear theory are shown in Figure 18, when Hj/Lo and S are com- 
bined into a single parameter. 
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Bese auaws cnoidal third-order Stokes 
cnoidal-linear 
——- linear 


Figure 12. Comparison of experimental and theoretical wave breaking 
angles (bottom slope, S = 0.10). 
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Table 2. Limit for validity of cnoidal 
theory for normal incidence. 


Bottom Slope S#0.0l 


Figure 13. Wave breaking angle, ap, as a function of 
deepwater wave angle, a,, and deepwater 
steepness, Hj/L>, on a bottom slope, S = 0.01. 
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Bottom Slope S$ = 0.02 


Figure 14. Wave breaking angle, ap, asa function of 
deepwater wave angle, a 5, and deepwater 
steepness, H,/Lj, on a bottom slope, S = 0:02. 


Bottom Slope S003 


Figure 15. Wave breaking angle, ap, as a function of 
deepwater wave angle, a,, and deepwater 
steepness, Ho/Lo> on a bottom slope, S = 0.03. 
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Bottom Slope S=0.05 


Figure 16. Wave breaking angle, ap, as a function of 
deepwater wave angle, a,, and deepwater 
steepness, H,/L5, on a bottom slope, S = 0.05. 


Bottom Slope S#O.1 


Figure 17. Wave breaking angle, op, as a function of 
deepwater wave angle, a,, and deepwater 
steepness, Ho/Lo> on a bottom slope, S = 0.10. 
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Figure 18. Prediction of wave breaking angles using 
linear theory. 


The remaining breaking wave characteristics are easily determined 
as shown in the following example: 
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(a) Cnoidal-linear 
From Figure 15 is found 
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Using Snell's law 


The wave height can be found from the breaking criterion 
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Finally, the depth of breaking is computed from the dispersion 
relationship, which is solved for 
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(b) Linear theory 
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IV. SUMMARY AND CONCLUSIONS 


It appears that no single wave theory accurately predicts the trans- 
formation of a wave from deep to shallow water. In shallow water, cnoidal 
theory successfully describes the shoaling of wave height but overestimates 
the wave celerity and wavelength. All previous studies emphasize an accurate 
determination of wave height, but little attention is paid to the wavelength. 


When considering longshore currents and sediment transport the angle 
of wave breaking becomes a parameter of high importance and, therefore, also 
the wave refraction process which is intimately connected with wavelength. 


In this study a "hybrid" wave theory is proposed, consisting of 
cnoidal wave height and linear wavelength transformation in the regime where 
cnoidal theory applies. In comparison with existing data this theory is 
found to predict wave height and wavelength better than previous theories, 
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especially for larger wave steepnesses. More definite conclusions must 
await new experimental data. This basic information seems to be missing 
in the existing literature. 


Nomographs are presented for determining the breaking wave 
characteristics for given deepwater characteristics and bottom slope for 
a plane beach. 


The remaining difficulties encountered in the shoaling and re- 
fraction of a monochromatic wave are further amplified when dealing with 
a more realistic directional spectrum of waves. What are the equivalent 
monochromatic breaking wave characteristics that produce the same sediment 
transport or longshore current? This question seems to be one of the more 
important to face in the future research. 
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APPENDIX 


Hyperbolic Wave Shoaling 


It is well known that near breaking surface gravity water waves are 
highly nonlinear and therefore not well characterized by linear theory. 
For prediction of breaking wave characteristics it is necessary to use 
higher order wave theories and, in particular, a higher order cnoidal wave 
solution could appear to be attractive because of the improved fit to data 
found when using first-order cnoidal theory. Iwagaki (1968) derived a 
practical asymptotic solution to second-order cnoidal waves called 
“hyperbolic waves." This approximation takes m = 1 and E(m) = 1 for 
K(m) > 3 where K(m) and E(m) are the complete elliptical integrals of the 
first and second kind. 


If the wave height transformation is written in the usual manner 
as a product between a shoaling and refraction coefficient 


H=H,K,-K (A-1) 


then Iwagaki derived 
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where 
d = water depth 
2 
L = a + na) = deepwater wavelength (third-order Stokes 
5 g wave) 
HS = unrefracted deepwater wave height 
K = complete elliptic integral 
H = wave height 
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and Iwagaki found by empirical curve fitting 
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Such that the elliptic integral, K, could be approximated by 
1/2 n m 
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—-2 9 fl-a@} (A-4) 
TVvg/d 
Finally, in hyperbolic waves the wavelength is given by 
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Again, it should be pointed out that equations (A-2) to (A-5) are only valid 
mone IC 2 Oe 


To evaluate the hyperbolic wave theory the shoaling factor K, 1s com- 
pared with data and first-order cnoidal theory (Brink-Kjaer and Jonsson, 
(1973). The results are shown in Figure A-1. 


It is evident that the results given by the hyperbolic theory are 
worse than those obtained using linear cnoidal theory. For large wave 
steepness the hyperbolic theory exhibits a decrease in K, near breaking 
similar to fifth-order Stokes waves which is also due to nonhomogeneous 
convergence of the perturbation series. Finally, the wavelength predicted 
by hyperbolic theory (not shown here) is quite close to the cnoidal wave- 
length which is in itself questionable as examined earlier. Further 
attempts at using hyperbolic wave theory to predict breaking wave character- 
istics were abandoned because of the above shortcomings. 
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